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ABSTRACT 

We investigate the Spitzer/IRAC properties of 36 z ~ 7 z 850 — dropout galaxies and 3 z ~ 8 Yogs 
galaxies derived from deep/wide-area WFC3/IR data of the Early Release Science, the ultradeep 
HUDF09, and wide-area NICMOS data. We fit stellar population synthesis models to the SEDs 
to derive mean redshifts, stellar masses, and ages. The z ~ 7 galaxies are best characterized by 
substantial ages (> 100 Myr) and M/Ly w 0.2. The main trend with decreasing luminosity is that 
of bluing of the far— UV slope from f3 ~ —2.0 to /3 ~ —3.0. This can be explained by decreasing 
metallicity, except for the lowest luminosity galaxies (0.1 L*_ 3 ), where low metallicity and smooth 
SFHs fail to match the blue far— UV and moderately red H — [3.6] color. Such colors may require 
episodic SFHs with short periods of activity and quiescence ("on-off" cycles) and/or a contribution 
from emission lines. The stellar mass of our sample of z ~ 7 star forming galaxies correlates with 
SFR according to logM* = 8.70(±0.09) + 1.06(±0.10) logSFR, implying star formation may have 
commenced at z > 10. No galaxies are found with SFRs much higher or lower than the past averaged 
SFR suggesting that the typical star formation timescales are probably a substantial fraction of the 
Hubble time. We report the first IRAC detection of Y 09S — dropout galaxies at z ~ 8. The average 
rest-frame U — V f» 0.3 (AB) of the 3 galaxies are similar to faint z ~ 7 galaxies, implying similar 
M/L. The stellar mass density to M vv ,ab < -18 is p*{z = 8) = l.S^o x 10 6 M Mpc" 3 , following 
logp*(z) = 10.6(±0.6) - 4.4(±0.7)log(l + z) [M s Mpc" 3 ] over 3 < z < 8. 
Subject headings: galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

Until recently, only a modest number of relatively 
bright z > 7 galaxies were known, mostly from wide- 
area NICMOS ^Bouwens et alJl200a lOesch et al.ll2009aL 
R. J. Bouwens et al. i n preparation) and gr ound- 
based searchesdOuchi et al.ll2009t iCastellano et alj|2009t 
iHickev et all I2009D . The arrival of WFC3/IR aboard 
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HST has dramatically improved the situation by iden- 
tifying large numbers of z > 7 galaxies by their red- 
shifted UV light. Here we report on z > 7 galaxies se- 
lected from the WFC3/IR Early Release Science (ERS) 
obser vations over the GOODS-South field ( Wilkins et al. 
120091 R. J. Bouwens et al. in preparation), comple- 
mented with candidates from the recent ul tradeep sur- 
vey with WFC3/IR over th e HUDF 09 field (I Oesch et al l 
1 20101 iBouwens et al.ll201~0"at se e also lMcLure et al.ll2009t 
Bunker et alJl2009UYan et alj|2009j ). 

Little is known about the stellar masses, metal pro- 
duction, and the contribution of star formation to reion- 
ization in these galaxies. Mid-infrar ed observations w ith 
the InfraRed Array Camera (IRAC: lFazio et al.ll2004[ ) on 
Spitzer can been used to constrain the stellar masses and 
ages, which has lead to the surprising discovery of quite 
massive ~ lQ 10 Mp, galaxies at z > 6 (|Evles et al J 120051 : 
lYan et aLI 120061 : IStark et aT1l2009D and app reciable ages 
(200-300Myr) and M/L s as early as z ~ 7 (lEgami et all 
MM . ILabbe et all [20061: IGonzalez et all Mm . The 
overall results suggest the galaxies formed substantial 
amounts of stars at even earlier times, well into the 
epoch of reionizat ion (jStark et alj|2007t lYan et alj|2006t 
ILabbe et al1l2010f ). 

In this Letter, we study the stellar populations of the 
largest sample of z > 7 galaxies with IRAC measure- 
ments to date, focusing on correlations with luminosity 
and stellar mass, and implications for the mass density 
and 2^8. We adopt an = 0.3, Sl\ — 0.7 cosmology 
with Hq = 70 km s~ 1 "Mpc~ 1 . Magnitud es are in the AB 
photometric system (|Oke fc Gundll983l ). 
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2. OBSERVATIONS AND STELLAR POPULATION 
MODELING 

Our sample consists of sources derived from the ultra- 
deep WFC3/IR HUDF, the deep WFC3/IR ERS, and 
wide-area NICMOS over the CDF-S and CDF-N. Can- 
didates were selected using the z ~ 7 zsro — dropouts an d 
z ~ 8 ln98— d r opouts , as u sed in IQesch et all 1 2010), 
IBouwens etaLl (|2010bD . and IGonzalez et al.l (|2009t ) fsee 
also R. J. Bouwens et al. in preparation). We now briefly 
discuss the IRAC photometry from the new ERS sample. 

The Spitzer/IRAC data over the WFC3/ERS area in 
the CDF-S (~ 23.3 hours integration time) was ob- 
tained from the Great Observatories Origins Deep Sur- 
vey (GOODS; M. Dickinson et al. in preparation). Ft . 
The IRAC depths in the 3.6 and 4.5— bands are 27.1 and 
26.5 magnitude (la, total, point source), respectively. 
Obtaining reliable IRAC fluxes of the candidates is chal- 
lenging because of the contamination from the extended 
PSF wings of nearby foreground sources. We remove con- 
taminating flux, by modeling the candidates and nearby 
sources using their isolated flux profiles and positions in 
the deep WFC3/IR maps as templates. We convolve the 
templates to match the IRAC PSF, simultaneously fit 
them to the IRAC map leaving only the fluxes as free 
parameters, and su btract the best- fit models to the fore- 
ground sources /seelLabbe et al.ll2006|:lWuyts et alJ l2007t 
IGonzalez et al.ll2009Hde Santis et al.H2007D . After clean- 
ing the IRAC images, we perform conventional aperture 
photometry in the 3.6 and 4.5 bands in 2'.'5 diameter 
apertures on 15 of the original 18 zsso - dropout galaxies 
over the WFC3/IR ERS. Three were too close to bright 
sources for reliable measurement. Fluxes were corrected 
by a factor xl.8 to account for light outside the aper- 
ture (consistent with point source profiles). Six of the 
15 galaxies are undetected in IRAC ([3.6] < 26.5,2c). 
Errors include the uncertainty in the best-fit confusion 
correction, added in quadrature. 

We complete the sample with 21 z ~ 7 galaxies with 
IRAC measuremen ts from the HUDF (jOesch et alJl2010t 
lLabbe et al.l I2010D and from t he re cent wide-area NIC- 
MOS search (jGonzalez et al.l 120091 R. J. Bouwens, in 
preparation). The total sample consists of 36 z ~ 7 
galaxies spanning 4 magnitudes in H 1(i Q. To investigate 
trends with magnitude, we stacked the flux densities of 
the galaxies in three ~ 1— magnitude bins centered on 
Hiso w 26, 27 and 28, containing 11, 15, and 10 galaxies 
respectively. The uncertainties are determined by boot- 
strapping. Stacking increases the SNR, in particular for 
faint galaxies, where the uncertainty in the mass is driven 
by the SNR in IRAC (see ?). 

We derive stellar masses and redshifts by fitting stellar 
populations synthesis models to the average SEP fluxes 
using the -fitting code FAST (jKriek et alJl2009Q . We 
adopt Bruzual fc Chariot! (|2003L BC03) models with a 
ISalpeter l (|1955D initial mass function (IMF) between 
0.1 — 100 Mq. We explore several SFHs and the effects 
of metallicity and dust. The diffe rences with more re- 
cent models (e.g., [Marastonl 120051 Chariot & Bruzual, 
in prepara tion) are small an d will not be considered in 
detail fsee lLabbe et al.lfeOlOt ). We fit models smoothed 
to a resolution of 100 A rest-frame, corresponding to 

11 This paper u ses data release DR2 of epoch 2, available from 
http: / /data. spitzer . caltech . edu/popular/goods/ 



the approximate width of the dropout selection win- 
dows. Adopting a lKroup a (2001) IMF reduces the stellar 
masses and SFRs by 0.2 dex, but does not change other 
parameters or the quality of fit. The typical uncertain- 
ties in the derived average stellar masses, SFRs, age, and 
Ay for the stacked SEDs are 0.15 dex, 0.25 dex, 0.3 dex, 
and 0.1 mag, respectively. The photometry and best- fit 
model parameters are presented in Table 1. 

3. STELLAR POPULATIONS AND STAR FORMATION 
HISTORIES AT Z ~ 7 

Figure 1 (left panel) shows the broadband SEDs of 
the z ~ 7 2^850 — dropout galaxies in the three magnitude 
bins, with the best-fit BC03 stellar population models. 
The overall SED shapes are remarkably similar, with a 
pronounced jump between -ffi6o and [3.6] (or rest-frame 
(U — V) « 0.5) indicative of a modest Balmer break and 
evolved stellar populations (> lOOMyr). Focusing on the 
far— UV continuum, we find the slope fx cx A* 9 (traced 
by the J125 — £fi60 color) to be very blue, and decreasing 

from 2 at H 160 ~ 26 t o ~ - 3.0 at H wo ~ 28. As 

discussed bv lBouwens et all (|2010bO . such extremely blue 
slopes require low dust content Ay < 0.1, very low metal- 
licities and/or very young ages. S mall 0.04 mag c hange s 
in the WFC3/IR zeropoints (e.g., iMcLure et all (2009)) 
would cause changes of f=s 0.17, comparable to the ran- 
dom uncertainties. The red -H160 ~ [3-6] color, however, 
implies more evolved stellar populations (> lOOMyr), 
leaving the models seemingly unable to match the en- 
tire SED. 

To explore the mismatch further we consider in more 
detail the effects on metallicity, SFH, and nebular emis- 
sion. Fig. 1 (right panel) shows the J125 — -Hi60 versus 
#160 — [3-6] colors of the observed stacked SEDs. The 
lines show predictions of 0.2 Z@ BC03 models for vari- 
ous continuous SFHs (rising, constant, declining). Gen- 
erally, evolved models are able to reproduce the joint 
J125 — Hieo versus Hiqq — [3.6] colors of the more lumi- 
nous z — 7 galaxies, but not the colors of the faintest, 
bluest galaxies. The arrows show the effect of changes in 
model assumptions, which we will discuss now: 

1) Metallicities: Low metallicities (e.g., 0.2 Z®) do 
a decent job of producing much bluer than Solar at 
a given -ffigo — [3.6] color. Very low metallicities (1/50 
Solar) produce even bluer 0's, but also bluen H 160 — 3.6. 
Metallicity alone appears not enough to fully resolve the 
discrepancy. 

2) Nebular Emission: Nebular emission lines (NEL) 
likely contribute to the IRAC fluxes, reddening the 
#160 — [3-6] color. Empirical estimates of [OIII]5007 
emission at z > 2 are scarce, but we can infer the 
possible effect from the observed strength of Ha at 
z ~ 2.2 (jErb et al.l 1200(1 . Assuming W Ha = 200 A and 
W [O ///]4959,5007+ff/3 = 2.5 x W Ha , ap p ropriate for 0.2 Z Q 



galax ies at z > 2 (|Erb et all l2006t iBrinchmann et all 
2008) , and adopting a redshift distribution of z = 6.9 ± 
0.5 (jOesch et al. 2010) we calculate a contribution of 0.18 
mag to [3.6] and 0.14 mag at [4.5] (see Fig 1). Note how- 
ever, that no measurements of nebular lines at z ~ 7 ex- 
ist. A second possible effect is nebular continuum emis- 
sion (NEC), which would cause a redde ning of both 
and Hieo — [3.6] (shown are the models of lSchaererll2002l 
with Z = 1/50 Z Q , t = 300 Myr CSF, and 0% escape 
fraction). Both would reduce the discrepancy between 
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Fig. 1. — (Top panels) Representative IRAC [3.6] image stamps of WFC3/ERS z ~ 7 2850 — dropouts (sorted by magnitude) after 
subtracting neighbors. Brighter sources may have lower SNR after accounting for confusion, (left panel) Stacked broadband SEDs from 
the combined NICMOS, WFC3/UDF and WFC3/ERS samples, averaged in ~ 1-mag bins centered on Hi 60 ~ 26,27 and 28. The data 
include HST/ACS, NICMOS, and WFC3/IR, groundbased K, and IRAC [3.6] and [4.5]. Upper limits are 2cr. The gray histogram shows 
the i?i6o-band magnitude distribution (the peak is 9 sources). The best-fit BC03 stellar population models at z = 6.9 are shown. The 
overall SED shapes are similar with a break between Hiqq and [3.6], expected f or evolved stellar popu lations (> lOOMyr). The far— UV 
slope (traced by J125 — #160 ) bluens towards fainter Higo magnitude (as found Bouwcns ct al. 2010b). (right panel) Comparison of the 
observed average J125 — Hiso versus -Hieo — [3.6] color (filled circles) with predictions of 0.2 Zq BC03 models for various SFHs (solid 
lines). Ellipsoids show la uncertainties. The SFHs are linearly increasing with time (blue), CSF (purple), episodic (light green) with a 
50% dutycyle and 40Myr duration, and exponentially declining (red) with r = lOO Myr. Also show n are the effects on the m odel colors 
of metallicity, reddening (Cal zetti et al. 2000), nebular continuum emission (NEC) (Schacrcr 2002), and line emission (NEL) IjErb et al.l 
2006; Brinchmann ct al. 2008). The episodic model is shown twice, once only with the effect of line emission (+ NEL) and once with both 
line and nebular continuum emission (+ NELC). The lines terminate at the 750Myr. 



models and data by allowing the models to extend to 
redder -£/i6o — [3.6] at a given j3. Note that the extremely 
blue j3 may require hi gh escape fractions f esc > 0.3 (e.g., 
iBouwense t al. 2010b), which would reduce the contribu- 
tion of nebular emission. 

3) Star Formation Histories: Declining SFHs can 
match the colors of the most luminous, redder galaxies, 
but their far— UV continua are too red for lower luminos- 
ity galaxies. In contrast, strongly rising SFRs (SFRoc t a , 
a > 1) exhibit blue /3 but never reach red i?i6o — [3-6] 
in a Hubble time (z = 7) and are formally excluded at 
95% confidence. CSF is an compromise, providing red 
H\m ~ [3-6] colors generated by on-going assembly of 
stellar mass, and blue far— UV continua from on-going 
star formation. Finally, episodic SFH with a 50% duty 
cyle and 40 Myr duration (i.e., 20 Myr "on", 20 Myr 
"off") are found to have an interesting mix of properties. 
The luminous active phase of the cycle produces a bluer 
far— UV continuum for a given -H160 — [3-6] color than 
CSF. The reverse is true in the dimmer passive state. 
The net result for a steep UV LF function is that the 



luminosity weighted average of cycling galaxies displays 
bluer i&i—UV at a given -H160 — 3.6 than CSF, also re- 
ducing the discrepancy. 

In summary, the model colors match the observations 
of luminous z ~ 7 galaxies reasonably except for the low- 
est luminosity galaxies, where the blue /3 ~ — 3 and the 
red -ffi6o — [3.6] ~ 0.6 colors remain challenging to fit. 
Low metallicity CSF models come close, but a contribu- 
tion from nebular line emission to the [3.6]— band and/or 
episodic SFHs are likely needed to resolve the mismatch. 

4. STAR FORMATION RATE VERSUS STELLAR MASS AT 

Z ~ 7 

Independent constraints on the SFHs can be obtained 
from the relation between SFR and stellar mass, as 
shown in Fig. 2. The SFRs are calculated from the 
monoc hromatic 1500A lum inosity following the prescrip- 
tion of iMadau et al.l ()1998l ) and corrected for dust using 
the best-fit Ay . The galaxies are grouped in bins of SFR 
(or M1500) centered on log SFR « 0.4, 0.8, 1.2. 

The stellar mass of our z ~ 7 sample correlates strongly 
with the SFR, producing logM* = 8.70(±0.09) + 
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Fig. 2. — Average UV— derived, dust-corrected SFR versus stel- 
lar mass for z ~ 7 galaxies. Black symbols denote galaxies from 
the NICMOS {circles), WFC3/IR ERS (squares) and WFC3/IR 
HUDF sample (trangles). The purple squares show the average 
stellar mass in bins of SFR centered on log SFR ?a 0.4,0.8,1.2. 
The purple triangles assume 0.2 and 0.15 mag contribution by 
emission lines to the [3.6] and [4.5] bands. The scatter in M* in 
bins of M1500 is ~ 0.3 dex. The diagonal lines show the maximum 
stellar mass CSF stellar population can form in a Hubble time 
(dotted) or since z = 10 (dashed). Galaxies with strongly increas- 
ing/declining SFRs would lie well above/below the lines, respec- 
tively, but few such systems are found. The inset shows simulated 
galaxies with random formation times and exponentially declining 
T = 200 Myr SFHs, showing a different distribution with larger 
scatter. (Top panels) Stacked images in the [3.6] band in bins of 
SFR. 

1.06(±0.10) log SFR and corresponding to a constant 
M/L v w 0.200 The derived stellar mass with emis- 
sion lines (see §3) would be lower by « 0.17dex. The 
scatter around the relation is fairly low 0.3 dex, but 
that does not exclude significant short-term SFR vari- 
ation, e.g., the episodic model ([J3]) predicts log M* = 
8.45 + f.Olog SFR 1500 and a scatter of ~ 0.3 dex. We 
find no galaxies with SFRs much lower or higher than 
the past averaged SFR (i.e., strongly bursting or sup- 
pressed. Such galaxies would have satisfied our dropout 
criteria and would lie in the upper left or lower right cor- 
ner in Fig. 2. Their absence suggests that that the typical 
star formation timescales are probably a substantial frac- 
tion of the Hubble time. Instead we find that only 4/22 
sources with log SFR > 0.7 are undetected at [3.6] and 
no galaxies in the sample have SFRs substantially less 
than M*/t HubMe . 
To illustrate the diagnostic power of the M* — SFR 

12 The M/L also be calculated from the i ndividual galaxies, 
yielding the same answer (Gonzalez ct al. 2009). 
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Fig. 3. — (top panels) Individual and stacked images of the 3 
2 ~ 8 Yo98 _ dropout galaxies in our WFC3/ERS sample, shown in 
inverted grayscale. Two of the 3 sources are individually detected 
in the IRAC [3.6]— band, (bottom panel) The average broadband 
SED of the 3 z ~ 8 galaxies and the best-fit stellar population 
model (solid blue line). For comparison, we show the best-fit model 
to z = 6.9 2850 dropouts of similar -Hieo magnitude (dashed green 
line), shifted by —0.5 mag. The overall shapes are similar, with 
the 2 ~ 8 galaxies being slightly bluer in -ffieo — [3.6] compared to 
2 ~ 7 galaxies. Upper limits are 2<r. 

diagram, we show in the inset in Fig. 2 a simulation of 
galaxies with random formation times and exponentially 
declining SFRs (r = 200 Myr) to the same selection lim- 
its as our observed sample. The distribution is clearly 
different, with no correlation between M* and SFR, sug- 
gesting that star formation timescales for z ~ 7 galaxies 
are probably longer than that. 

5. THE STELLAR MASS DENSITY AT Z ~ 8 

The detection of z ~ 8 galaxies with IRAC is enticing 
as it enables us to place stronger constraints on the stel- 
lar masses of the highest redshift galaxies than possible 
from the iax—UV alone. Recent studies in the HUDF 
have found no dete ction (individual o r stacked) for z ~ 8 
galaxy candidates (|Labbe et al. 2010), leaving estimates 
of the stellar mass density at these redshifts highly un- 
certain. 

Here we perform photometry on the 3 z ~ 8 
!o98 — dropout galaxies in the WFC3/ERS sample of R. J. 
Bouwens et al. (in preparation). These candidates are 
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Fig. 4. — The evolution of the integrated stellar mass den- 
sity. The red circle shows the z m 7.7 mass density, derived 
from the integrated UV— luminosity density of (Bouwens et al., 
in preparation) and the mean M/L derived here. All data are con- 
verted to a common limit Mjjv,ab < —18, using the UV LFs of 
IBouwens etaLl H2010al) . The WFC3/ERS data are corrected by 
+0.42 dex. The z = 3 - 7 luminous samples (M UViAB < -20) 
from the literature ([Stark ct al. 2009 blue circles, Gonzal ez et al.l 
(2009), magenta circle) are corrected by +0.38, +0.46, +0.57, and 
+0.75 dex at z =4, 5, 6, and 7, re spectively. The HUDF requires 
no correction l ILabbe et alj 12010. green circle). T he red square 
shows the z 8.2 sample of (Bouwens ct al. 2010a), but using the 
M/L derived here for the z « 7.7 sources. The dashed line shows 
a cx (1 + z)~ 4,4 evolution. The floating error bar indicates the 
expected cosmic variance for the z ~ 8 sample. 

brighter than found in the HUDF {H « 27 versus H « 
28). Two galaxies are detected at [3.60 and we calculate 
an average SNR = 3.7 in [3.6] and SNR = 1.6 in [4.5] for 
the stack of all three (see Fig 3, top panels). The best fits 
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mass-to-light ratios M/Ly = 0.15 and M/L1500 = 0.1, 
and log SSFR — —8.7. Overall these properties are com- 
parable to H sa 27 z ~ 7 Zs50 — dropouts, suggesting mod- 
est evolution in the M/L bet ween z = 8 and z = 7 . 

Following t he ap proach of iGonzalez et al.l (|2009D and 
ILabbe et all (l2010h . we derive integrated stellar mass 
densities at z ~ 8 by multiplying the UV— luminosity 
densit ies integrated to Muv,ab = — 18 ([Bouwens et al.1 
l2010al R .J .Bouwens et al. in preparation) by the 
mean M/L derived for the z ~ 7.7 galaxies, yielding 
p*(z = 8) = 1.8±i 7 a x 10 6 M Q Mpc" 3 . We also re- 
compu te the z ~ 8 stellar mass density of ILabbe et al.l 
(2010) using the same M/L. Figure 4 shows the evolu- 
tion of the stellar mass density from z = 3 to z = 8. 
The evolution over 3 < z < 8 is well approximated by 
logp*(z) = 10.6(±0.6)-4.4(±0.7)log(l + z) [M Q Mpc" 3 ] 
over 3 < z < 8. 

13 We caution that source ERSy-2376440061 at [3.6] is close to 
pixels that are affected by "Muxbleed" , which we subtracted using 
a 3rd order polynomial fit to the 20" X 20" background before 
performing photometry. 

Following Labbe ct al. (2006) we report SFH weighted age I, . , 
where age m = t/2 for CSF and t is the time elapsed since the start 
of star formation. 



6. SUMMARY 

Using a large sample of 36 2s5Q — dropout galaxies based 
on deep/wide-area WFC3/IR data from the Early Re- 
lease Science, ultradeep data from the HUDF09, and 
wide-area NICMOS programs, we investigate the stellar 
population properties at extreme redshifts z > 7. The 
main results are: 

• The average rest-frame far— UV slope at z ~ 7 
becomes bluer with decreasing luminosity, from 
P 2.0 fi:_,) to ~ -3.0 ( 0-1 £* =3 ), as re- 
ported bv IBouwens et al.l (|2010bl ). The rest- frame 
U — V becomes bluer as well, but is still moder- 
ately red U — V s» 0.3 at 0.1 L* z=3 , apparently 
excluding extremely young ages < lOOMyr. If the 
ages inferred from the simple model fits are cor- 
rect, galaxies started forming stars very early-on, 
perhaps as high as z > 10. The blue far— UV 
slope and red U — V colors remain a challenge to 
fit, however, even for sub-solar metallicity models. 
Episodic SFHs with periods of activity and quies- 
cence and/or a (~ 0.2 mag) contribution of emis- 
sion lines to the [3.6]— band may be required to 
resolve the mismatch. 

• The derived stellar masses correlate with the SFRs 
at z ~ 7 according to logM* = 8.70(±0.09) + 
1.06(±0.10)log SFR, with relatively low scatter 
~ 0.25 dex. Emission line contributions of ps 0.2 
mag to both [3.6] and [4.5] would shift the re- 
lation by s» —0.2 dex in mass. The absence of 
galaxies with SFRs much lower or higher than the 
past averaged SFR (i.e., strongly bursting or sup- 
pressed) suggests that that the typical star for- 
mation timescales are probably a substantial frac- 
tion of the Hubble time. Note that instantaneously 
quenched galaxies may fade too quickly to be se- 
lected as dropout galaxies. 

• The first Spitzer/IRAC detection of z = 8 galax- 
ies and their red average i?i6o — [3.6] ~ 0.55 sug- 
gest that luminous early galaxies may have sub- 
stantial M/Ly k. 0.15, similar to z — 7 galax- 
ies. The derived stellar mass density then in- 
creases gradually with time following logp*(z) = 
10.6(±0.6)-4.4(±0.7)log(l + z) [M Mpc" 3 ] over 
3 < z < 8. 

Deeper IRAC data on z > 7 galaxies are needed to 
bolster these results. More nebular emission line mea- 
surements of z > 2 galaxies would help to understand 
the possible contribution to the broadband fluxes. Ad- 
ditional modeling of the distribution of SFR versus M* 
is needed to decipher the SFHs of z ~ 7 galaxies. Larger 
samples would enable a more secure assessment of the 
mass density evolution beyond z = 8. 
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Star formation rates and stellar masses to z = 7 — 8 
TABLE 1 

Summary of photometry and modeling of z ~ 7 — 8 dropout galaxies 
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SEDs of z ~ 8 Yogs-dropouts 
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Average colors of 2 ~ 7 2850-dropouts 
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SFRs and stellar masses of 2 ~ 7 2gso-dropouts 
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Note. — The optical— to— near-IR fluxes arc measured in 0"4 diameter apertures. Spitzcr/IRAC fluxes are measured on the confusion- 
corrected maps in 2'/ 5 diameter apertures. Fluxes are corrected to total assuming point source profiles. Units are nanojy for the SEDs 
and AB magnitudes for the average colors. The total sample consists of 15 new sources from the WFC3/ERS sample (Bouwens et al., in 
preparation), 9 sources from the NICMOS sample (UDF-1417,964,GNS-1,2,3,4,5,CDFS-4627,HDFN-1216; Gonzalez ot al. 2010), and 12 
from the WFC3/UDF sample (UDFz-4471, 4257, 3955, 3958,3722,4314,3677,3744,4056,3638,3973,3853; Oosch ct al. 2010, Labbc ct al. 2010). 
The stacked Y — band of the z ~ 7 galaxies is a c ombination of the Y ips and Yo98 bands. SFR1500 is the SFR. derived from the 1500A 
monochromatic luminosity using the prescription of Madau ct al. (1998) and corrected for dust using the best-fit Ay (mean < Av >— 0.13 
mag). The stellar masses arc based on BC03 0.2Zq models using exponentially declining SFHs with 8 < log r < 11 and reddening A\r < 0.3. 
Uncertainties are determined by bootstrapping. 



